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Abstract 

A  mechanistic  model  that  interprets  the  oxidation  behavior  of  the  diborides  of  Zr,  Hf  and  Ti  in  the  temperature  range  of  ^1000- 
1800  °C  was  formulated.  Available  thermodynamic  data  and  literature  data  for  vapor  pressures  and  diffusivities  were  used  to  evaluate 
the  model.  Good  correspondence  was  obtained  between  theory  and  experiments  for  weight  gain,  recession  and  scale  thickness  as  func¬ 
tions  of  temperature  and  oxygen  partial  pressure.  At  temperatures  below  about  1400  °C,  the  rate-limiting  step  is  the  diffusion  of  dissolved 
oxygen  through  a  film  of  liquid  boria  in  capillaries  at  the  base  of  the  oxidation  product.  At  higher  temperatures,  the  boria  is  lost  by 
evaporation,  and  the  oxidation  rate  is  limited  by  Knudsen  diffusion  of  molecular  oxygen  through  the  capillaries  between  nearly  columnar 
blocks  of  the  oxide,  MO2. 
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1.  Introduction 

Components  in  high-flow  environments,  such  as  the 
leading  edge  of  a  hypersonic  vehicle,  are  subjected  to  very 
high  temperatures  and  require  that  the  geometric  integrity 
be  maintained  during  service  [1,2].  The  high  temperatures 
and  the  low-pressure  environment  make  oxidation  and 
evaporative  loss  critical  factors  in  material  selection  [3]. 
The  diborides  of  Zr  and  Hf,  which  are  termed  ultra-high- 
temperature  ceramics  (UHTCs),  have  high  melting  points 
and  are  among  the  most  promising  materials  for  long-term 
service  under  hypersonic  conditions  [3,4].  The  diborides  are 
also  attractive  due  to  their  high  thermal  conductivity, 
which  is  useful  to  offset  the  stringent  thermal  gradients 
imposed  by  aerothermal  heating.  These  diborides  have  also 
been  considered  as  coatings  for  protection  of  carbon-based 
composites  [5]  and  titanium  alloys  [6]. 


Corresponding  author. 

E-mail  address:  Triplicane.parthasarathy@wpafb.afmil  (T.A.  Partha¬ 
sarathy). 


A  history  of  studies  on  the  oxidation  of  diborides  is  given 
by  Opeka  et  al.  [3].  Based  on  the  promise  shown  from  early 
works,  a  good  oxidation  screening  study  of  the  refractory 
diborides  of  Ti,  Zr,  Hf,  Nb  and  Ta  was  conducted  by  Kauf¬ 
man  et  al.  [7,8]  and  later  summarized  by  Fenter  [9].  They 
reported  that  the  diborides  of  Zr  and  Hf  are  the  most  oxida¬ 
tion  resistant.  They  further  reported  that  additions  of  SiC 
up  to  20%  by  volume  further  improved  the  oxidation  resis¬ 
tance  of  these  materials.  Opeka  et  al.  [3,4,10]  pointed  out 
that  the  high  melting  points  and  low  vapor  pressures  of 
the  oxides  and  sub-oxides  of  Zr  and  Hf,  and  the  beneficial 
vapor  pressure  vs.  Pq^  relationships  for  BO^  species,  are 
responsible  for  the  superior  high-temperature  resistance  of 
these  materials.  Using  volatility  diagrams,  they  further 
showed  that  the  diborides  of  Zr  and  Hf  appear  to  have  the 
least  problem  with  respect  to  disruption  of  the  product  scale 
by  vapor  species.  Zr  and  Hf  carbides  are  limited  by  CO  par¬ 
tial  pressures  that  exceed  1  atm  above  1730  °C,  yielding  a 
porous  and  non-pro tective  oxide  scale.  Silica  has  very  low 
diffusivities  for  oxygen  and  is  a  very  protective  scale  at  more 
moderate  temperatures,  but  active  oxidation  to  gaseous  SiO 
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disrupts  the  scale  above  approximately  1800  °C  [3,4,10]. 
Consequently,  SiC,  Si3N4  and  transition  metal  silicides  are 
not  very  useful  at  ultra-high  temperatures.  In  contrast,  boric 
oxide  is  glassy,  flows  to  fill  pores  in  the  oxide  scale  and 
exhibits  a  1  atm  partial  pressure  (for  gaseous  B2O3)  only 
at  1950  °C,  despite  its  well-known  tendency  to  evaporate 
at  moderate  temperatures  (above  ^1200  °C).  The  gaseous 
B2O3  evaporates  at  open  surface  sites  without  disrupting 
the  scale.  When  SiC  is  added  to  the  system,  further  improve¬ 
ments  are  realized  since  the  silica-based  glass  offers  greater 
resistance  to  evaporation  than  B2O3  alone  [11]. 

Thus  far,  the  understanding  of  the  superior  oxidation 
behavior  of  the  diborides  has  been  primarily  qualitative. 
For  engineering  design  and  applications,  a  quantitative 
model  is  needed  that  predicts  all  aspects  of  oxidation.  It 
is  desirable  to  predict  factors  such  as  scale  thickness,  sub¬ 
strate  recession  and  weight  change,  under  a  set  of  complex 
conditions  of  temperature  and  environment.  Many  experi¬ 
mental  works  have  been  conducted  under  isothermal  con¬ 
ditions  and  in  a  laboratory  environment,  while  real 
engineering  conditions  involve  thermal  gradients,  environ¬ 
mental  gradients  and  turbulent  gas  flow. 

The  long-term  objective  of  this  work  is  to  develop  a 
model  that  can  be  applied  to  aid  engineering  design  for 
hypersonic  environments.  However,  the  complex  nature 
of  the  materials  makes  it  necessary  to  start  with  a  model 
for  simpler  systems  and  conditions,  where  reliable  and 
well-defined  experimental  data  are  available.  In  this  paper, 
we  present  a  model  that  is  found  to  predict  the  scale  thick¬ 
ness,  metal  recession  and  weight  changes  for  the  oxidation 
of  Zr,  Hf  and  Ti  diborides  under  isothermal  and  slow  gas 
flow  conditions.  Various  factors  such  as  residual  boria  con¬ 
tent  and  dependencies  on  temperature  and  oxygen  partial 
pressure  are  also  predicted  and  found  to  be  in  reasonable 
agreement  with  experiments.  This  model  provides  a  reason¬ 
able  basis  for  future  work. 


2.  Conceptual  framework 

The  model  reported  in  this  work  is  based  on  concepts 
of  the  oxidation  process  derived  from  micro  structural 
details  available  in  the  literature.  The  morphology  of 
the  oxidation  product,  as  well  as  the  phases  present  and 
their  distribution,  has  been  reported  by  several  researchers 
[3,4,9,10,12].  For  ZrB2  (or  HfB2),  oxidation  in  air  results 
in  a  dense  adherent  oxide  scale  consisting  of  two  phases, 
Zr02  (or  Hf02)  and  B2O3.  No  intermediate  phases  or 
borates  of  these  systems  are  observed  at  the  interfaces. 
The  distribution  of  the  phases  in  the  scale  varies  with  tem¬ 
perature.  At  lower  temperatures  (<1000°C),  a  glassy 
B2O3  film  is  observed  on  top  of  the  (Zr02  +  B2O3)  scale, 
but  this  external  boria  layer  is  absent  at  higher  tempera¬ 
tures.  At  all  temperatures  a  porous  zirconia  scale  is  one 
component  of  the  oxidation  product.  At  low  and  interme¬ 
diate  temperatures,  the  pores  in  the  zirconia  are  filled  with 
boria.  Further,  the  micro  structure  of  the  Zr02  appears  to 
change  from  equiaxed  grains  at  temperatures  below 
~1000  °C,  to  columnar  grains  at  higher  temperatures. 
At  higher  temperatures,  the  outer  surface  of  the  scale  is 
nodular  Zr02,  and  the  cross-section  reveals  a  columnar 
structure  of  Zr02  with  a  glassy  B2O3  in  between  the 
columnar  grains.  A  similar  transition  in  oxidation  rates 
has  been  noted  with  respect  to  oxygen  partial  pressure. 
The  dependence  of  the  oxidation  rate  on  the  oxygen  par¬ 
tial  pressure  changes  from  linear  to  no  dependence  as  the 
temperature  is  increased  above  1150  °C  for  ZrB2  [13].  In 
addition  to  these  two  experimentally  observed  tempera¬ 
ture  regimes,  it  is  expected  that  a  high-temperature  regime 
exists  where  there  is  very  little  or  no  boria  in  the  scale. 
The  three  regimes  are  shown  schematically  in  Fig.  1.  In 
this  work,  we  focus  on  the  intermediate  temperature 
regime,  which  is  the  regime  of  engineering  interest.  The 
model  developed  here  predicts  the  temperature  range  of 


Intermediate  Temp  (-1000  to  -1800  °C) 

«>.aaXa 

:  ZrB, 

AA.AAA 

-  $00^ 

^2^3(9) 

ZrOj' 

Fig.  1.  The  oxidation  products  formed  during  oxidation  of  ZrB2  in  three  temperature  regimes. 
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the  intermediate  regime  and  the  oxidation  behavior  within 
this  intermediate  regime. 

The  conceptual  framework  for  the  oxidation  in  the 
intermediate  temperature  regime  is  shown  schematically 
in  Fig.  2a,  and  the  steps  considered  in  the  model  are  shown 
in  Fig.  2b.  For  oxidation  at  intermediate  temperatures,  oxi¬ 
dation  of  the  metal  (in  the  diboride)  results  in  an  evapora¬ 
tion-resistant  refractory  oxide  that  forms  a  porous 
skeleton.  The  oxidation  of  the  boron  results  in  a  glassy 
boric  oxide  that  flows  to  fill  the  base  of  the  porous  skele¬ 
ton.  At  the  surface,  the  boria  evaporates.  Both  oxidation 
reactions  occur  at  the  refractory  substrate-scale  interface 
by  the  inward  diffusion  of  oxidant  species.  At  steady  state, 
oxidation  takes  place  by  both  gaseous  diffusion  of  oxygen 
through  the  open  pores  in  the  oxide  skeleton,  and  by  dis¬ 


solved  oxygen  diffusing  through  the  liquid  boria  (at  the 
base  of  the  porous  skeleton  beneath  the  surface)  to  reach 
the  ZrB2-scale  interface.  The  two  reactions  are: 

ZrB2  “h  2^^  — ^  Zr02  B2O3  (1) 

B2O3  (1)  ^  B2O3  (g)  (2) 

The  rate  of  dissolution  of  oxygen  into  B2O3  is  not  taken 
to  be  rate  limiting  in  the  model  developed  below.  The 
transport  of  oxygen  through  the  refractory  oxide  is 
assumed  to  be  negligible  compared  to  that  through  the 
pores  filled  with  glassy  boria.  At  temperatures  of  interest 
(above  1000  K),  B2O3  (v)  has  been  shown  to  be  the  domi¬ 
nant  vapor  species  at  oxygen  partial  pressures  greater  than 
that  in  equilibrium  with  ZrB2  and  Zr02  [12]. 


a 


Scale  thickness,  L 


Fig.  2.  (a)  Schematic  sketch  of  mechanisms  involved  oxidation  of  ZrB2  in  air,  in  the  intermediate  temperature  regime  (1000-1800  °C).  (b)  Schematic  of  the 
mechanistic  steps  considered  in  the  model. 
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3.  The  model 


The  objective  of  the  model  is  to  predict  the  kinetics  of 
scale  thickness,  substrate  recession  and  weight  change  as 
a  function  of  time  at  temperature  and  oxygen  partial  pres¬ 
sure.  A  schematic  for  the  model  is  shown  in  Fig.  2b.  Dense 
Zr02  is  assumed  to  be  a  barrier  to  oxygen  permeation,  due 
to  its  negligible  rate  of  ambipolar  diffusion.  The  validity  of 
this  assumption  is  presented  in  the  Appendix.  The  porosity 
between  Zr02  columns  is  assumed  to  be  continuous,  and 
tortuosity  is  neglected. 

The  ambient  air  is  taken  to  consist  of  N2  and  O2  only. 
A  fraction,  /,  of  the  scale  is  taken  to  be  porous  and  pro¬ 
vide  a  continuous  pathway  for  gaseous  diffusion  and  for 
the  evaporation  of  liquid  B2O3  at  the  surface,  which  is 
taken  to  be  a  perfect  sink  for  B2O3.  The  relative  rates 
of  filling  the  pore  by  oxidation  at  the  substrate-scale 
interface  and  the  evaporation  of  B2O3  determines  the 
extent  to  which  the  pores  are  filled  with  B2O3.  At  a  given 
time,  t,  the  total  scale  has  a  thickness  of  L  and  its  inner 
portion  is  filled  with  liquid  B2O3  to  a  thickness,  h.  The 
distance  L-h  is  the  diffusion  distance  in  the  pore  over 
which  gaseous  oxygen  must  diffuse  to  reach  the  B2O3. 
The  oxygen  is  assumed  to  dissolve  in  B2O3  and  continue 
to  diffuse  to  the  substrate-scale  interface  where  the  oxi¬ 
dation  reaction  occurs.  As  explained  in  the  Appendix, 
the  dense  zirconia  columns  are  considered  to  be  impervi¬ 
ous  to  oxygen. 

At  the  substrate-scale  interface,  reaction  (1)  takes  place. 
Thus,  at  steady  state,  the  diffusional  flux  of  oxygen  must  be 
balanced  by  the  formation  rates  of  Zr02  and  B2O3. 

\J02\  =  2^^6203  =  2^Zr02-  (3) 


The  flux  of  gaseous  oxygen  in  the  porous  channel  that  is  L- 
h  long  is  given  by: 


\Jo,\=fDo, 


^02  ^02 

L-h 


(4) 


where  D  is  the  diffusion  coefficient,  C  the  concentration,  the 
superscript  “i”  refers  to  the  interface  between  B2O3  (1)  and 
the  gas  phase,  which  is  at  a  distance  h  from  the  substrate- 
scale  interface,  and  the  superscript  “a”  refers  to  the  scale- 
ambient  interface.  Similarly,  the  flux  of  B2O3  (v)  from  the 
surface  of  B2O3  (1)  to  the  ambient,  is  given  by: 


‘/B2O3  — /^B203 


^B203  ^B203 

L-h 


(5) 


Taking  JB2O3  =  ^B203  (quasi-steady  state),  and  combining 
Eqs.  (3)-(5)  the  concentration  of  molecular  oxygen  at  the 
B2O3  liquid-vapor  interface  is  obtained  as: 

^02  ~  ^02  “  2  LAq  ^  ^^®203  “  ^6203)-  (^) 

The  concentration  at  the  ambient  surface  is  taken  to 
be  zero.  The  concentration  Cg203  B2O3  (v)  at  this 

interface,  which  does  not  depend  on  the  local  oxygen  pres¬ 
sure,  is  obtained  from  vapor  pressure  data  [14]: 


Cb,03  =  (atm)  =  3  x  10*  exp  . 

(7) 

The  permeability  of  oxygen  through  the  liquid  B2O3 
depends  on  the  oxygen  activity  gradient  across  the  liquid 
layer  of  thickness  t.  From  the  ideal  gas  law,  the  oxygen  par¬ 
tial  pressure  at  the  B2O3  liquid-vapor  interface  is  given  by. 


Po,=RTC'o,- 


(8) 


The  activity  of  oxygen  at  the  scale-B203  interface  is  given 
by: 


K 


reaction  1 


<^ZrB2  (-^02)^'^^ 
<^B203<^Zr02 


(9) 


Using  the  thermodynamic  data  of  Barin  [14]  for  reaction 
(1),  and  equating  the  activities  of  the  solid  phases  to  unity, 
Eq.  (9)  gives  the  following  for  the  oxygen  activity,  Fq^,  at 
the  substrate  ZrB2-B203  interface  (s): 


(atm)  =  5  X  10^^  exp 


99,967\ 
T  )' 


(10) 


The  values  given  by  Eq.  (10)  are  consistent  with  the 
maps  presented  by  Fahrenholtz  [12].  The  diffusivity  of 
oxygen  through  boria  has  been  measured  by  Tokuda 
et  al.  [15],  who  found  a  linear  dependence  of  the  perme¬ 
ation  flux  on  oxygen  partial  pressure  and  deduced  that 
molecular  oxygen  diffuses.  The  oxygen  flux  across  the  li¬ 
quid  boria,  ^02(^203),  is  related  to  the  partial  pressure 
gradient  through  the  oxygen  permeability  coefficient,  77, 
as  [16,17]: 


1^02(6203)1 


7702-B2O3 

h 


(11) 


The  oxygen  permeability  coefficient,  77o2  B2O35  is  obtained 
from  the  literature  [15,17]  as: 


7702-B2O3  (mol/m  s  atm)  =  0.15  exp  .  (12) 

Substituting  Eq.  (12)  into  Eq.  (11)  and  equating  the  flux  of 
oxygen  through  the  gaseous  layer,  given  by  Eq.  (4),  an 
equation  relating  the  boria  layer  thickness,  A,  to  the  scale 
thickness,  L,  is  obtained: 


_  _ fl?^02  B2O3  ) _ 

^  “  Do,{Cl^  -  +  7702-0203(7^^  -  ' 

(13) 

The  equations  for  the  rate  of  change  of  scale  thickness  (L) 
and  mass  ( W)  of  the  scale  are  given  by  accounting  for  the 
flux  balance  using  Eq.  (3) 


dL 

Tt 


1  \  (  1  \  dlTzr02 
1  ~  f )  VPZr02/  *1^ 


1  \1  Mzr02 
1-/75  Pzr02 


(14) 
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where  p  refers  to  density  and  Mto  molar  volume.  Combin¬ 
ing  with  Eq.  (4),  one  obtains: 


dL  2/Mz.oA/  f\  Cq,  -  Cq, 
~APzrojV-f)  L{\-q)- 


(15) 


Integration  of  Eq.  (15)  gives  a  parabolic  equation  for 
growth  of  the  scale: 


=  2 


2  fM zx02\ 
^  \  Pzr02  / 


^02 

(1  -Q) 


t. 


(16) 


The  magnitude  of  recession  of  the  substrate  is  given  by: 


^Zr02/Pzr02 

The  total  weight  change  per  unit  area  is  given  by: 


(17) 


AW 


^Pzr02  ~  f)  ^  PB2O3 


^PzrB2  • 


(18) 


It  remains  to  evaluate  the  diffusion  coefficient  of  gaseous 
oxygen  through  the  porous  region  filled  with  oxygen, 
B2O3  (g)  and  nitrogen.  The  diffusion  coefficient  in  a  mul¬ 
ti-component  gaseous  system  can  be  approximated  by 
[18]: 


1 

EfA/AA’ 


(19) 


where  D  refers  to  diffusivity,  n  to  mole  fraction,  subscript 
“1,  (2,  . . .,  /)”  refers  to  the  diffusivity  of  species  1  in  a  med¬ 
ium  of  i  species,  subscript  “1  —  /”  refers  to  the  diffusivity  of 
species  1  in  a  binary  mixture  of  species  1  and  i.  The  diffu¬ 
sivity  Di_2  of  species  1  in  a  binary  gas  mixture  with  species 
2  is  given  by  [19]: 


D\2 


0.0018583r3''V(lM)  +  O/M2) 

r\2QT>P 


where  D\_2  is  the  gas  diffusivity  in  cmVs; 


(20) 


^  _  1.06036  0.193000  1.03587 

D  -  ro.15610  +exp(o.47635r)  ‘^exp(1.52996r) 

1.76474 

'^exp(3.8941ir)’ 

T’*  =  ,  ei2  =  ri2  =  0.5(ri  +  rj)  in  A,  P  =  pressure 

in  atm;  M,  =  molecular  weight  (g/mol). 

The  parameters  needed  for  the  above  expressions  can  be 
obtained  from  the  work  of  Svehla  [20].  When  the  size  of  the 
pore  is  small,  the  mean  free  path  is  longer  than  the  charac¬ 
teristic  dimension  of  the  system,  and  the  diffusion  of  gases 
is  governed  by  Knudsen  diffusion  [21]: 


A 


4  /87?r\  r 

3V^y  2’ 


(21) 


Z)eff  =  (77k‘+77r-2)''-  (22) 

Using  Eqs.  (19)-(22),  Eqs.  (16)-(18)  give  the  scale  thick¬ 
ness,  recession  and  weight  gain  as  a  function  of  tempera¬ 
ture,  time  and  partial  pressure  of  oxygen. 

4.  Model  predictions  compared  with  literature  data 

4.1.  ZrBz 

The  model  was  verified  by  comparing  its  predictions  with 
experimental  data  reported  in  the  literature  for  the  oxidation 
of  ZrB2.  As  can  be  seen  from  the  development  in  Section  3, 
most  of  the  parameters  needed  for  the  prediction  can  be 
obtained  from  the  literature.  Two  parameters  in  the  model 
that  are  not  well  known  are  the  pore  fraction  and  pore 
radius,  and  they  are  likely  to  vary  between  experiments 
depending  on  starting  microstructure  and  oxidation  condi¬ 
tions  such  as  temperature  and  environment.  The  values  of 
0.05  and  0.5  pm,  respectively,  were  chosen  to  get  the  best 
fit  with  all  of  the  experimental  data  for  ZrB2,  but  were  kept 
the  same  for  HfB2  and  TiB2.  From  what  can  be  inferred  from 
micrographs  of  oxide  scales  presented  in  the  literature,  these 
values  are  quite  plausible.  The  sensitivity  of  the  predictions 
to  the  choice  of  these  variables  was  examined  and  is  pre¬ 
sented  in  a  later  section.  The  other  significant  unknowns 
are  the  velocities  and  water  contents  of  the  atmospheres  in 
the  experiments.  Future  experimenters  are  encouraged  to 
report  these  parameters;  they  can  have  large  effects. 

Extensive  thermogravimetric  analysis  (TGA)  data  on 
weight  change  vs.  time  in  250Torr  of  pure  oxygen  was 
reported  by  Tripp  and  Graham  [22].  They  compared  their 
results  with  the  parabolic  rate  constant  data  from  Berko- 


0.00055  0.00075  0.00095 

1  /  Temperature  (K'^) 

Fig.  3.  The  experimental  data  from  three  different  works  in  the  literature 
[12,22,23]  on  the  parabolic  rate  constant  of  measured  weight  gain  of 
samples  in  250  Torr  of  pure  oxygen  at  different  temperatures  are  shown 
compared  with  the  model  predictions.  The  porosity  in  the  scale,  /,  was 
taken  to  be  0.05  and  the  pore  radius  to  be  0.5  pm. 


where  is  the  Knudsen  diffusivity,  M  the  molecular  mass 
of  the  diffusing  species  and  r  the  radius  of  the  porous  path¬ 
way.  The  effective  diffusivity  is  given  by  [21]: 
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witz-Mattuck  [13]  and  Kuriakose  and  Margrave  [23].  Fig.  3 
shows  the  results  from  all  three  experimental  works,  com¬ 
pared  with  the  model  predictions,  for  an  assumed  pore 
fraction  of  0.05  and  pore  radius  of  0.5  pm.  The  weight  gain 
data  of  Tripp  and  Graham  are  lower  than  the  model  pre¬ 
dictions  at  temperatures  above  1200  °C.  The  lower  values 
of  experimental  data  could  arise  from  a  low  gas  flow  rate 
used  in  the  experiment,  while  the  model  assumes  a  perfect 
sink  for  B2O3  (v).  However,  when  the  model  is  compared 
with  all  of  the  reported  data,  the  correspondence  is  taken 
to  be  good  in  that  it  matches  well  at  low  temperatures 
and  is  a  good  representation  of  the  widely  scattered  data 
at  high  temperatures. 

Fenter  summarized  the  metal  recession  of  several  dibo¬ 
rides  in  air,  as  a  function  of  temperature  [9].  The  results 
for  ZrB2  in  Fig.  4a  are  compared  to  the  model  prediction. 
The  model  compares  very  well  with  data  up  to  ^1850  °C, 
above  which  the  model  predicts  that  the  high-temperature 
regime  dominates  and  that  all  of  the  boria  will  evaporate 


as  soon  as  it  forms.  The  experimental  data  show  a  signifi¬ 
cant  enhancement  in  recession  at  this  temperature.  The 
model  predicts  the  transition  temperature  correctly.  The 
reason  for  the  upward  trend  in  recession  is  not  clear;  one 
possibility  is  significant  spallation  of  MO2  due  to  the 
absence  of  boria  at  higher  temperature. 

Finally,  Opeka  et  al.  [3,10]  and  Talmy  et  al.  [24]  have 
reported  the  mass  change  and  oxide  layer  thicknesses  of 
ZrB2  oxidized  in  air  or  Ar/02  mixtures  with  an  oxygen  par¬ 
tial  pressure  of  0.2.  Their  data  are  compared  with  the 
model  predictions  in  Fig.  4b  and  c.  Again,  the  correspon¬ 
dence  is  taken  to  be  quite  reasonable  given  the  assumed 
geometrical  factors. 

4.2.  HfB2 

The  model  is  easily  extended  to  other  diborides  by  sim¬ 
ply  using  the  appropriate  thermodynamic  data  and  the 
appropriate  physical  properties  such  as  molecular  weight. 


1800 °C  1600 °C 


1  /  Temperature  (K"'*) 


1  /  Temperature  (K'^)  1  /  Temperature  (K'^) 

Fig.  4.  (a)  Model  predictions  of  the  recession  in  air  as  a  function  of  temperature  is  shown  compared  with  data  of  Kaufman  and  Clougherty  [7,8]  and 
Fenter  [9].  (b)  The  model  predictions  for  mass  change  are  shown  compared  to  experimental  data  obtained  in  Ar/02  gas  in  2  h  [10].  (c)  The  scale  thickness 
obtained  after  5  h  in  air  by  Opeka  et  al.  is  shown  compared  to  the  predictions.  The  porosity  in  the  scale,  f,  was  taken  to  be  0.05  and  the  pore  radius  to  be 
0.5  pm. 
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Fig.  5.  Comparison  of  model  predictions  with  experimental  data  on  oxidation  of  HfB2.  (a)  Recession  in  flowing  air  from  Kaufman  and  Clougherty  [7,8] 
and  Fenter  [9].  (b)  Parabolic  rate  constant  for  weight  gain  as  a  function  of  temperature  from  Berkowitz-Mattuck  in  250  Torr  of  oxygen  [13].  The  porosity 
in  the  scale,/,  was  taken  to  be  0.05  and  pore  radius  to  be  0.5  pm. 


density,  etc.  For  HfB2,  one  obtains  the  oxygen  partial  pres¬ 
sure,  at  the  HfB2-B203  interface  as: 

P^^(HfB2)  (atm)  =  6  x  10^  exp  .  (23) 

Two  reports  in  the  literature  have  characterized  the  oxi¬ 
dation  behavior  of  HfB2.  Fenter  [9]  summarized  the  mea¬ 
sured  recession  rates  in  flowing  air,  while  Berkowitz- 
Mattuck  [13]  reported  on  the  parabolic  rate  constants  in 
250  Torr  of  pure  oxygen.  These  data  are  compared  with 
the  model  predictions  using  a  pore  fraction  of  0.05  and 
pore  radius  of  0.5  pm  in  Fig.  5.  The  predicted  recession  is 
higher  than  the  measured  values,  while  the  predicted  para¬ 
bolic  rate  constants  are  lower. 

4.3.  TiB2 

For  TiB2,  one  obtains  the  oxygen  activity,  Fq^,  at  the 
TiB2-B203  interface  as: 

P^^(TiB2)  (atm)  =  5  x  10'°  exp  (^-  .  (24) 

By  use  of  the  appropriate  physical  properties  for  Ti02  and 
TiB2  the  oxidation  behavior  was  predicted,  once  again 
assuming  a  pore  fraction  of  0.05  and  a  pore  radius  of 
0.5  pm.  The  predictions  were  evaluated  by  comparing  with 
experimental  data  available  on  weight  change  and  reces¬ 
sion  reported  by  Koh  et  al.  [25].  This  comparison  is  shown 
in  Fig.  6. 

5.  Parametric/sensitivity  studies 

A  parametric  study  was  conducted  to  assess  model  sen¬ 
sitivity  and  to  determine  the  effects  of  experimental  vari¬ 
ables.  Fig.  7  shows  the  temperature  dependence  of  the 
parabolic  rate  constant  for  oxidation  of  ZrB2.  A  clear 


Fig.  6.  Comparison  of  model  predictions  with  experimental  data  on  TiB2. 
(a)  Weight  gain  with  time  in  air  and  (b)  scale  thickness  after  oxidation  in 
air.  The  porosity  in  the  scale,/  was  taken  to  be  0.05,  with  a  pore  radius  of 
0.5  pm.  The  data  were  obtained  from  Koh  et  al.  [25]. 
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1400  °c 


Fig.  7.  Results  of  a  parameteric  study  on  the  effect  of  temperature  on  oxidation  behavior  of  ZrB2.  (a)  The  parabolic  rate  constant  and  (b)  the  relative 
fraction  of  liquid  boria  in  the  scale.  The  model  predicts  a  change  in  the  apparent  activation  energy  for  the  parabolic  rate  constant  that  correlates  with  loss 
of  boria  from  the  scale. 
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Fig.  8.  The  effect  of  pore  radius  (Knudsen  effect)  on:  (a)  recession,  (b)  boria  fraction  in  scale,  (c)  oxide  thickness  and  (d)  parabolic  rate  constant. 
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change  in  apparent  activation  energy  is  seen  at  around 
1400  °C.  The  figure  includes  the  predicted  fraction  of  boria 
(/z/L)  in  the  scale.  The  change  in  activation  energy  clearly 
arises  from  the  loss  of  boria  from  the  scale. 

Next,  the  model  was  used  to  study  the  effect  of  varying 
the  pore  fraction  and  pore  radius,  to  assess  its  sensitivity. 
These  parameters  are  difficult  to  measure  and  have  not  been 
reported  experimentally.  Further,  since  the  model  neglects 
the  tortuosity  of  the  porous  scale,  the  pore  radius  and  pore 
fraction  used  in  the  model  must  be  considered  to  be  “effec¬ 
tive”  parameters  rather  than  the  actual  physical  values.  It  is 
quite  possible  that  these  parameters  themselves  vary  with 
temperature  or  other  experimental  conditions  such  as  envi¬ 
ronment.  Thus  it  is  useful  to  conduct  a  sensitivity  study. 

Fig.  8  shows  the  effect  of  pore  radius  on  the  recession, 
scale  thickness,  parabolic  rate  constant,  and  the  thickness 
fraction  of  boria  in  the  scale.  The  choice  of  pore  radius  is 
significant  when  the  pore  radius  is  smaller  than  about 
1  pm.  However,  the  effect  is  insignificant  at  temperatures 
below  about  1400  °C,  which  corresponds  to  the  tempera¬ 
ture  at  which  boria  evaporation  becomes  significant. 


Fig.  9  shows  the  effect  of  pore  fraction  on  the  oxidation 
behavior  of  ZrB2  in  air.  The  rate  constant,  scale  thickness 
and  recession  all  increase  by  about  an  order  of  magnitude 
when  the  pore  fraction  is  increased  from  0.025  to  0.2.  Thus 
the  effect  is  nearly  linear. 

Comparing  Figs.  8  and  9  to  Figs.  3-5  it  can  be  seen  that 
very  close  correspondence  between  the  model  and  experi¬ 
mental  data  could  be  obtained  if  pore  fraction  and/or  pore 
radius  are  taken  to  vary  rather  than  remain  constant  for  all 
scales.  The  scatter  in  the  experimental  results  alone  implies 
that  these  geometrical  parameters  may  vary  depending  on 
the  exposure  temperature,  and/or  experimental  conditions 
such  as  flow  rates  of  gases,  and  humidity,  and  they  may 
even  change  somewhat  during  the  course  of  an  experiment. 
The  literature  data  is  typically  incomplete  with  respect  to 
details  of  the  environment,  especially  humidity,  and  the 
scale  micro  structures. 

Fig.  10a  summarizes  the  effect  of  oxygen  partial  pressure 
in  O2/N2  gas  mixture  of  1  atm  total  pressure,  on  the  para¬ 
bolic  rate  constant  as  a  function  of  temperature.  The 
dependence  is  linear  at  1000  °C  and  becomes  negligible  at 


1  /  Temperature  (K'^) 


Fig.  9.  The  effect  of  pore  fraction  on:  (a)  rate  constant,  (b)  oxide  thickness  and  (c)  recession,  for  a  change  in  pore  fraction  from  0.025  to  0.2. 
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Oxygen  partial  pressure,  atm 


Oxygen  partial  pressure,  atm 
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Fig.  10.  Parametric  studies  of  the  dependence  of:  (a)  parabolic  rate  constant  on  the  oxygen  partial  pressure  in  O2/N2  gas  mixture  (1  atm  total  pressure). 
The  P02  dependence  itself  is  dependent  on  (b)  the  pore  size  and  (c)  pore  fraction. 


1600  °C;  this  is  consistent  with  the  experimental  observa¬ 
tions  of  Berkowiz-Mattuck  [13].  Fig.  10b  shows  that  the 
pore  radius  is  predicted  to  have  a  significant  effect  at  low 
oxygen  partial  pressures  and  the  pore  fraction  has  a  near- 
linear  effect  (Fig.  10c)  as  might  be  expected. 

6.  Summary 

A  physical  model  to  predict  the  oxidation  behavior  of 
refractory  diborides  in  the  intermediate  temperature  regime 
(1000-1800  °C)  is  presented.  The  model  assumes  that  the 
refractory  oxide  product  does  not  support  significant  oxy¬ 
gen  diffusion  due  to  a  low  oxygen  vacancy  concentration 
and  to  a  lack  of  sufficient  electronic  conductivity.  The 
porosity  in  the  base  of  the  oxide  is  assumed  to  be  filled  with 
liquid  boria  which,  however,  evaporates  from  the  surface. 
A  comparison  of  the  model  with  experimental  data  in  the 
literature  shows  reasonable  agreement  in  predicting  weight 
change,  metal  recession,  oxide  scale  thickness  and  the  tem¬ 
perature  dependence  of  the  parabolic  rate  constant  for 
ZrB2,  even  with  postulated  values  for  pore  geometry  and 
unknown  experimental  flow  rates.  Comparison  with  lim¬ 
ited  literature  data  for  HfB2  and  TiB2  show  that  the  model 
may  have  a  general  applicability  to  refractory  diborides. 


although  data  on  oxygen  diffusion  and  electronic  conduc¬ 
tivity  of  the  oxides,  especially  at  high  temperatures,  will 
be  required  to  refine  the  model. 
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Appendix.  Discounting  ambipolar  oxygen  diffusion  in 
zirconia 

In  the  model  for  refractory  diboride  oxidation,  the  diffu¬ 
sion  of  oxygen,  via  ambipolar  diffusion,  is  considered  to  be 
negligible  because  of  the  low  magnitudes  of  oxygen  vacan¬ 
cies  and  p-type  electronic  conduction  in  the  zirconia  oxida¬ 
tion  product.  In  the  relatively  high  oxygen  activities 
encountered  in  the  oxidation,  only  p-type  electronic  con- 
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duction  needs  to  be  considered;  n-type  conduction  is  negli¬ 
gible.  In  fact,  data  are  not  available  for  the  magnitude  of 
the  p-type  electronic  conductivity  in  relatively  pure 
Zr02,  but  some  measurements  have  been  made  for  calcia- 
stabilized-zirconia,  e.g.  Zro.85Cao.15O1  35,  as  a  function  of 
temperature.  The  interpretation  of  defect  equilibria  by  Las¬ 
ker  and  Rapp  [26]  permits  a  general  evaluation  of  p-type 
conductivity  as  a  function  of  dopant  concentration  and 
oxygen  activity  to  provide  magnitudes  for  the  ionic  and 
p-type  electronic  conductivities  for  compositions  with 
much  lower  dopant  concentrations. 

Consider  the  defect  equilibrium  for  zirconia: 

0^  +  2h^  to2  (g)+Vo,  (Al) 

where  Oq  is  an  oxygen  ion  on  an  oxygen  lattice  site,  Vq  is  a 
doubly  ionized  oxygen  vacancy  and  h  is  a  positive  hole. 
From  the  law  of  mass  action: 


Km 


[Vo]^ 

[h]^ 


(A2) 


If  the  zirconia  were  doped  by  substitution  of  an  aliovalent 
solute,  e.g.  Ca^^  ions,  on  the  Zr"^^  lattice  sites,  then  the  sim¬ 
plified  electrical  neutrality  condition  would  read: 

2[Ca"J=2[V6]«[h].  (A3) 

Upon  substitution  of  Eq.  (A3)  into  Eq.  (A2),  two  predic¬ 
tions  arise: 


[Vq]  (X  [Ca^j.]  7^f{Po2)  (A4a) 

[h]  ^  (A4b) 

From  the  usual  justified  assumption  of  concentration- 
independent  mobilities  of  the  oxide  ions  and  positive  holes, 
the  ionic  conductivity  and  p-type  electronic  conductivity 
obey  the  same  dependencies  as  Eqs.  (A4a)  and  (A4b). 

Patterson  [27]  has  presented  an  evaluation  of  the  ionic 
transference  number  for  Zro.85Cao.15O1. 35  (CSZ)  as  a  func¬ 
tion  of  temperature  to  establish  the  electrolytic  domain. 
From  a  knowledge  of  the  ionic  conductivity,  provided  by 
Patterson  et  al.  [28]  and  Patterson’s  transference  numbers 
[27],  the  p-type  electronic  conductivity  for  CSZ  can  be  cal¬ 
culated.  Then,  using  Eqs.  (A4a)  and  (A4b),  the  functions 
can  be  applied  to  nearly  pure  Zr02,  e.g.  containing 
150  ppm  divalent  impurity  concentration. 

From  Patterson  et  al.  [28]  for  0.15CaO-stabihzed  Zr02 
(CSZ),  using  the  ionic  conductivity  at  1273  K  and  the 
reported  activation  energy,  the  oxygen  partial  pressure- 
independent  ionic  conductivity  is  given  as: 

(Tion  (Q“‘  cm-')  =  2323 exp  •  (A5) 


Similarly,  the  p-type  electronic  conductivity  in  CSZ  can  be 
derived.  The  reported  hole  conductivity  is  2  x 
IQ-"'  Q-'  cm-'  at  1273  K  and  Po^  =  10*^  with  =  0.99. 

Further,  the  ionic  conductivity  equals  the  hole  conductivity 
at  (1/r,  P02)  values  of  (0.0003,  10^)  and  (0.0012,  10^"^). 


From  these  data,  and  combining  with  Eq.  (A5),  for  condi¬ 
tions  where  dion  =  the  hole  conductivity  is  obtained  as: 

(Q-'  cm-' )  =  23, 280^'’  exp  •  (A6) 

From  Lasker  and  Rapp  [26],  the  dependence  on  [Ca]  is  gi¬ 
ven  as: 


(7^  oc  [Cay'''^  (Tion  oc  [Ca"j'.  (A7) 

Combining  Eqs.  A5,  A6,  A7  one  obtains: 


(Tion  (i2  ‘cm  ')  =^^5^^2323 exp 


123, 500 Y 
8.314r  y’ 


(A8) 


dj:^  (Q  ^  cm  ^) 


/^[Cdopant]\ 


1/2 

23,280Po/^exp 


235,200\ 
8.314r  )' 
(A9) 


Consider  a  Zr02  composition  doped  with  150  ppm 
of  a  divalent  impurity.  At  1273  K,  the  ionic  and  hole 
conductivities  are  2  x  10“^  cm“^  and  1.64  x  10“^ 
cm“\  respectively,  with  a  transference  number  ^ion  = 
0.9918.  Thus  the  assumption  of  negligible  ambipolar  diffu¬ 
sion  is  validated  for  this  temperature. 

At  1873  K,  the  ionic  and  hole  conductivities  are 
8.35  X  10“"^  cm“^  and  2.03  x  10“"^  cm“\  respec¬ 
tively,  with  a  transference  number  ^ion  =  0.804.  Thus,  at 
these  higher  temperatures,  p-type  conductivity  becomes 
more  important  and  ambipolar  diffusion  of  oxide  ions  can¬ 
not  be  neglected.  However,  at  high  temperatures,  the  pores 
in  the  zirconia  are  expected  to  be  devoid  of  the  boria  glass 
due  to  evaporation.  Thus  it  is  unlikely  that  ambipolar  oxy¬ 
gen  diffusion  could  contribute  a  flux  that  is  competitive  with 
respect  to  gaseous  diffusion  of  oxygen  molecules  in  the  gas 
phase.  The  assumption  of  at  least  150  ppm  of  aliovalent 
impurity  in  the  zirconia  oxidation  product  is  not  unrealistic. 
The  remaining  issue  with  this  analysis  is  the  assumption  that 
the  lightly  doped  zirconia  product  would  exist  in  a  CaF2-type 
lattice  structure  for  the  scale;  however,  if  the  zirconia  scale 
existed  on  a  different  crystalline  state  (e.g.  tetragonal  or 
monoclinic),  the  oxygen  diffusivity  would  certainly  be  negli¬ 
gible  regardless  of  the  electronic  conductivity. 
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